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Abstract 

An algebraic model of Boson realization is proposed to study the vibrational 
spectra of a tetrahedral molecule, where ten sets of boson creation and annihilation 
operators are used to construct the Hamiltonian with symmetry. There are 
two schemes in our model. The first scheme provides an eight-parameter fit to 
the published experimental vibrational eigenvalues of methane with a root-mean- 
square deviation 11.61 cm~^. The second scheme, where the bending oscillators 
are assumed to be harmonic and the interactions between the bending vibrations 
are neglected, provides a five-parameter fit with a root-mean-square deviation 12.42 
cm~^. 
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I. INTRODUCTION 



The characterization of highly excited vibrational states has become one of the central 
goals in chemical physics. There were two general methods used to describe molecular 
vibrations. In the traditional approach the molecular Hamiltonian was parametrized in 
terms of internal coordinates [1] . The potential was modeled by the force field constants 
with many parameters due to poor knowledge of the large number of force constants. 
Those parameters have to be determined by fitting the spectroscopic data phenomeno- 
logically. 

As an alternative, an algebraic approach has been proposed for the study of polymolec- 
ular spectra. The first step toward the establishment of an algebraic approach was given 
by lachello, Levine and their co-workers [2] with the vibron model, where the rotation- 
vibration spectra of diatomic molecules are described in terms of a ^(4) algebra. Although 
this model was extended [3] to polyatomic molecules by introducing a -^(4) algebra for each 
bond, it is rather difficult to apply when the number of atoms in the molecule becomes 
larger than four [4] . 

Recently, an alternative technique [5] for the automatic computation of symmetrized 
local mode basis functions was used to provide a four parameter potential model for the 
stretching modes of octahedron XYq that gave an excellent fit to the published experi- 
mental vibrational eigenvalues of SFq, WFq and UFq. An improved algebraic method [6], 
where the one-dimensional Morse oscillator was described by the Lie algebra u{2), was 
proposed to provide another better fit to those experimental data with four parameters 
plus one fixed parameter that describes the anharmonic property of the Morse potential. 
This algebraic approach was extended to study the vibrational spectra, both stretching 
and bending, of tetrahedral molecules [4], where the interactions between stretching and 
bending vibrations were neglected and seven adjustable parameters plus two fixed param- 
eters Aj were used to fit the experimental data. This algebraic method was also used to 
study the vibrational spectroscopy and intramolecular relaxation of benzene [7] . 

In this paper we propose another algebraic model, the boson realization model, to 
study the vibrational spectra of tetrahedral molecules, where 10 coupled onc-dimcnsional 
anharmonic oscillators are described by 10 sets of bosonic creation and annihilation op- 
erators. The interbond interactions and the interactions between stretching and bending 
vibrations are expressed by the invariant combinations of the products of one creation 
operator and one annihilation operator such that the total number of vibrational quanta is 
conservative. The symmetrized bases arc used to simplify the calculation. There are two 
schemes in our model. The first scheme provides an eight-parameter fit to the published 
experimental vibrational eigenvalues of methane better than the previous results. The re- 
sults show that the bending oscillators are near harmonic ones, the interbond interactions 
between bending vibrations are quite weak, and the interactions between stretching and 
bending vibrations arc strong. Those conclusions reflect the properties of the molecular 
structure of methane. From the properties we propose our second scheme where the bend- 
ing oscillators are harmonic and there is no interaction between the bending vibrations. 
The second scheme provides a five-parameter fit to the experimental data of methane with 
the root-mean-square deviation 12.42 cm~^. It may be a model with the least parameters 
that well fits the published experimental vibrational data of methane. 

To some extent, our method is a generalization of that used in Ref. [5]. 
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This paper is organized as follows. In Sec. II the vibrational Hamiltonian of a tetra- 
hedral molecule is introduced in terms of ten sets of bosonic operators. In Sec. Ill the 
vibrational functions are combined into the symmetrized bases belonging to given rows 
of given irreducible representations of T^. In these symmetrized bases the Hamiltonian 
becomes a block matrix with eight parameters. The spurious states are ruled out in the 
calculation. In Sec. IV, fitting the published 19 experimental data [4, 8] for methane with 
the total number of quanta f < 3, in our first scheme with eight parameters we obtain 
the root-mean-square deviation of energy to be 11.61 cm~^. Our second scheme provides 
a good five-parameter fit. For comparison, the experimental data, the previous calculated 
energies by the algebraic model [4] and the present calculated results for methane ofv < 3 
in two schemes are listed in Table 2. The remaining calculated energies for methane up to 
V — 3 hy this boson realization model in the first scheme are also presented. The higher 
energy levels can be calculated straightforwardly. In Sec. V we give some conclusions. 



II. HAMILTONIAN 



We begin with enumerating 10 oscillators for an XY4 tetrahcdral molecule such as that 
in Fig. 1. The atom X is located at the center O of the tetrahedron, and four atoms Y at 
its vertices A, B, C, and D. The coordinate axes x, y, and z point from O to the centers 
of edges AC, AD, and AB, respectively. The first four equivalent oscillators describe the 
fundamental stretching modes {Ai © F2), and the other six equivalent ones describe the 
fundamental bending modes. As is well known [9], there are only five degrees of freedom 
for the bending vibrations {E © F2), so that the six bending oscillators must contain a 
spurious one. The spurious states related to the spurious degree of freedom should be 
ruled out in the later calculation. 
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Fig.l Schematic representation of a XY4 tetrahedral molecule. 



Now, for the tenoscillators we introduce 10 sets of bosonic operators and a^, 1 < 
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a <10, that satisfy the well known relations 

ap] = [al , al] = 0, 

(2.1) 



[tta , = 6al3, [aa , = , O^] = 0, 

\n) = \ni,n2, - ■ ■ ,nw), 



al\n) = y/Ua +1 I • • • , ria-l, {Ua + 1), ria+l, ,, 

aa\n) = ^/n^\■■■,na-l,{na-l),na+l,■■■), 

where \n) denotes the common eigenstate of the phonon number operators with the 
eigenvalues ria, respectively. 

Na = aiaa, \n) = n„ \n). (2.2) 

The first four bosonic operators a] (or a^), 1 < J < 4, describing the stretching 
vibrations, are the irreducible tensor operators belonging to the representations Ai © F2 
of Td- The other six a^^ (or a^), 5 < < 10, describing the bending vibrations, are those 
belonging to the representations Ai ® E ® F2. 

The energy of each oscillator depends upon the phonon number. For simplicity we 
assume that all oscillators are the Morse ones with two parameters uj and x, so that the 
energy of the ath oscillator can be expressed in the operator form: 

Es{N,) = N, {us - Xs (iV, + 1)} , 1 < J < 4 

Eb{N^) = {cUb - Xb{N^ + l)}, 5 < /X < 10 ^ ""^^ 

where the subscript s denotes the stretching vibration and b the bending one. The null 
energy has been removed. Although the Morse potential is known [10] to be not very 
suitable for all anharmonic oscillators, the deviation can be described by some more 
parameters that become important for higher energy levels. It was pointed out by lachello 
and Oss [7] that the Poschl- Teller potential is much more appropriate than the Morse one 
to describe the bending vibrations. However, the expression for the eigenvalues of the 
bound states is identical for two potentials [7]. 

As usual, neglecting the mixture of the states with different total number of phonons, 
and assuming to take the interactions up to the second order, we can express the inter- 
active potentials as the combinations of the products of one creation operator and one 
annihilation operator. The character table, the representation matrices of the generators, 
and the Clebsch-Gordan coefficients of the group were explicitly given in Ref. [4]. From 
that knowledge, there are obviously only five invariant combinations in addition to the 
sum of phonon number operators. The Hamiltonian now can be expressed in terms of the 
bosonic operators as follows: 

4 10 4 10 

7 { ^ 

+ A3 ^ (aj,a^+3 + 0)^+30^) + A4 < al ^ (a^ - a^+3) 

/i=5 I /i=5 

/ 8 \ ' 

+ aU as - ^ + og + aio I + ^ («2m - 02^-1) 

+ a\ ^-ag - ae + ^ a^ - aioj + //.C.j + A5 | [y.^ a] j a^ | + H.C. 

(2.4) 
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It will be seen in the next section that the term with A5 relates only to the spurious states 
so that it is not interesting to us. Removing this term, we obtain the Hamiltonian con- 
taining eight parameters. The term with A4 describes the interactions between stretching 
and bending vibrations. In the previous algebraic approach [4] ten parameters were in- 
troduced to fit the experimental data, where two spectroscopic constants A^i and N2, that 
are equal to cUs/xg — 1 and cub/xb — 1 in our notations, were taken as fixed parameters, one 
constraint was assumed to reduce one parameter, and the interactions between stretching 
and bending vibrations were neglected. In this meaning their algebraic approach [4], in 
comparison with our model, introduced three more parameters (-B12, -65,6, and -85,10), 
fixed two parameters A^i and N2, reduced one parameter by a constraint, and neglected 
one parameter A4 describing the interaction. 



III. IRREDUCIBLE BASES 

Since the Hamiltonian has symmetry, each eigenfunction can be combined such that 

it belongs to a given row of an irreducible representation of T^. The states that belong 
to the same irreducible space as partners must correspond to the same energy [11]. This 
degeneracy of the partners is called normal one. A symmetric perturbation never splits 
a normal degeneracy. The calculation for energy levels will be greatly simplified if those 
bases are used. 

The vibrational state of a tetrahedral molecule is described by the phonon numbers 
Ua of the ten oscillators. The first four numbers nj describe the stretching vibrations, and 
the next six numbers describe the bending vibrations. Those states can be combined 
into the irreducible bases belonging to given rows of given irreducible representations, 
respectively. For the general vibrations, those states with both stretching and bending 
vibrations should be further combined. 

Now, we discuss the combinations of the states for pure stretching vibrations. Briefly 
denote 1 711712^3^4000000) by \abcd), where the vanishing (/x > 5) are neglected in this 
notation. Firstly assume that a, b, c and d are all different from each other. Under the 
transformations of there are 24 independent states that span the regular representation 
space of Td- Through the standard group theory method [11] they can be combined into 
orthogonal bases belonging to ten irreducible representations: (f){Ai,abcd), (f){A2,abcd), 
(p^J\E, abed), <pf\E, abed), ^^(i?^, abed), (p^^\F2, abed), (p'^\F2, abed), abed), (t>f\F^, abed), 

and 0{Fi,abed). 

Similarly, for pure bending vibrations, briefly denote [00 007^57^6n77^87^97^lo) by \abcdef), 
where the vanishing rij {j < 4) are neglected in this notation. Two kinds of states should 
not be confused: one with four numbers describes the stretching vibration, and the other 
with six numbers describes the bending vibration. When a, b, e, d, e and / are all different 
from each other, we also have 24 independent states, spanning the regular representa- 
tion space of Td- Combine them into the irreducible orthogonal bases, and denote them 
by '4>{Ai,abedef), il){A2,abedef), ^pl^\E, abedef), ijj^^E, abcdef) , il>^\F2,abedef) , 
^(,2)(F2, abedef), ^l^\F2, abedef), tpi^^Fi, abedef), i)^^\Fi, abcdef), and ipi^^Fi, abedef). 
The explicit combinations of 0i/(r, abed) and ipiy{T, abedef) can be obtained from us upon 
request. 
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For example, for the fundamental stretching vibrations {v — 1) the irreducible bases 
are listed as follows: 



0(^,1000) = {|1000) + |0100) + |0010) + |0001)} 
0S'^(F2,1OOO) = 2-1 {jlOOO) - lOlOO) + lOOlO) - |0001)} 

0^^^(^2,1000) = {|iooo) - |oioo) - |oioo) + loool)} ^ ■ ^ 

0^^^(^2,1000) = 2-1 {|1000) + |0100) - |0010) - |0001)} 

Similarly, the irreducible bases of the fundamental bending vibrations {v = 1) are: 

V^(A, 100000) = 6-1/2 {1 100000) ^ loiOOOO) + looiooo) + lOOOlOO) 

+ lOOOOlO) + lOOOOOl)} 
^f^E, 100000) = (2\/3)-i {2|100000) - |010000) - |001000) 

+ 2|000100) - lOOOOlO) - lOOOOOl)} 
V'?^(£;, 100000) = 2-1 {1010000} - |001000) + lOOOOlO) - lOOOOOl)} ^^-^^ 
^'^^^2, 100000) = 2-1/2 iioiooOO) - 1000010)} 
^/;^^^(F2, 100000) = 2-1/2 ||ooiOOO) - |000001)} 
V'i^^(F2, 100000) = 2-1/2 iiioooOO) - |000100)} 

In those bases the Hamiltonian H given in (2.4) becomes a block matrix with the 
submatrices H(r,v), where F denotes the irreducible representation, and v is the total 
phonon number. Obviously, H{E, 1) is a 1 x 1 submatrix, but H{Ai, 1) and H{F2, 1) are 
2x2 submatrices: 

H{E,l) = a;6 - 2x6 - A2 + A3 

H(A ^) ^ ( '^s-2a;^ + 3Ai 2y/6\5 \ 

^ ) - [ 2VQX5 - 2x5 + 5A2 + A3 j (3.3) 

fT(p 1^ ^ f ujs-2xs-\i 2^2X4 \ 

^ ^ V 2x/2A4 uj, - 2x, - A2 - A3 j 

Note that the state ip{Ai, 100000) represents the fundamental spurious state and should 
be ruled out. 

In the traditional approach, the higher excited states are calculated by symmetrizing 
the products of the fundamental vibrational states [9]. However, it may be more easy 
to use the irreducible bases, subtracting the spurious states, for calculating the excited 
states. In the following, we calculate the excited states for the case with total phonon 
number f = 2 in detail. It is straightforward to calculate higher excited states in this way. 

Before calculation, we have to study an important problem of how to remove the 
spurious states. In the recent papers we see two methods of removal. lachello and Oss 
[7] placed the spurious states at the energies > 10 times the energies of the physical 
states by the projection operators. This method of removal is exact for harmonic bending 
vibrations and acquires a small error for anharmonic one. Instead, Lemus and Frank [4] 
directly eliminated the spurious states from both the space and the Hamiltonian. They 
demanded the matrix elements of the Hamiltonian related with the fundamental spurious 
state ip{Ai, 100000) vanishing. In our notation, they introduced a constraint in addition 
to A5 = that was assumed in Ref.[4] at beginning: 

uJb-2xb + 5X2 + X3 = (3.4) 
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so that H{Ai, 1) in (3.3) contains only one non- vanishing elements. 

In the present paper we develop the second method of removal. First of all, it seems 
to us that the constraint (3.4) is not necessary and reasonable because it restricts only 
the energy of the fundamental spurious state to be vanishing, but the energies of all 
other spurious states are non- vanishing. Secondly, we have to answer the problem how to 
identify the spurious states. Generally speaking, a state is identified as a spurious state if 
it contains ^'(^ii 100000) as a factor. In Ref. [4] (p. 8327) the simply additive definition 
for the product of two functions is used: 

\n) \m) = \{n + m)) 

\{n + m)) = |(ni + mi), (ria + m2), • • • , (nio + mio)) 

When two states \n) and \m) describe a pure stretching vibration and a bending vibration, 
respectively, (3.5) is correct. However, in the calculations of removing the spurious states, 
both states describe bending vibrations where (3.5) may not be suitable. 

Borrowing the idea from Ref. [7] , we want to find an identification rule for the spurious 
state such that the spurious species are separated, if possible, from the physical species 
in the matrix form of the Hamiltonian. As is well known, in the formalism of the boson 
reahzation the states \n) contains a factor (n!)~^/^: 

\n) = (n!)-i/2 i^^tj" |o) 
Therefore, we embed a factor in the definition (3.5) for product: 

WH- n{^^^f (3.6) 

\{n + m)) = I (ni + mi), (n2 + ma), •• -,(7110 + 77110)) 

In terms of this definition we calculate the general spurious states and find that the 
off-diagonal elements of the Hamiltonian between the spurious species and the physical 
species linearly depend upon Xb and A5, namely, under the conditions that Xf, = and 
A5 = the spurious species are totally separated from the physical species in the matrix 
form of the Hamiltonian. 

The condition A5 = is acceptable because it only appears in those off-diagonal 
elements. The condition x^ = means that the bending vibrations are harmonic. Fortu- 
nately, to our knowledge, in the known results Xf, is quite small {e.g. see Ref. [6], [12] and 
our results below). It is interesting to notice that the first method [7] of removal is exact 
also only for harmonic bending vibrations. 

Now, we turn back to calculate the excited states with v = 2. When v = 2, the stretch- 
ing vibrational states are separated into 5 sets: 0(Ai,2OOO), (I)^^\F2,2000), 0(Ai, 1100), 
(f)l}\E , 1100) , and 0^^^(^2,1100), and the bending vibrational states are separated into 
10 sets: ^(Ai, 200000), ^«(E, 200000), 7/^^(^2,200000), ^(Ai, 100100), ^W(^> 100100), 
i/j^^^Ai, 110000), il^l^^E, 110000), 7/;W(F2, 110000), 7/»(^HF2, 110000), and 7/;W(Fi, 110000). 
For the mixture states ^,y(r e Fi F2,f) of stretching and bending vibrations with 
t; = 2, we have to combine the stretching states (f),y{ri,v — 1) and the bending states 
'0j,(r2,w = 1) by the Clebsch-Gordan coefficients of T^i. 

In terms of the definition (3.6), direct calculation shows that there are 5 sets of spurious 
states with v — 2, belonging to the following irreducible representations: two Ai, one 
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E, and two F2. In the following we list only one spurious state for each irreducible 
representation space: 

0(^,1000) ^(Ai, 100000), 

200000) +-0(^1, 100100) + 2V'(Ai, 110000)} /V3 = 100000) }^ 

{'il^i^\E, 200000) + '02^^(£;, 100100) + 'ilji^\E, 110000)} /Vs 

= '0(A,iooooo)V'?^(£;, 100000), 

4^^(F2, 1000)V'(^, 100000), 

{il^i^\F2, 200000) + V2ipi^\F2, 110000)} /V3 = -0(^1, 100000)V'i^^(F2, 100000). 

(3.7) 

Removing the spurious states, we obtain the physical states belonging to given irre- 
ducible representations. There are five states belonging to representation Ai, five sets of 
states belonging to E, seven sets belonging to F2, and three sets belonging to Fi. Only 
one state for each irreducible representation space is listed in the following: 

/i(A,2) = 0(Ai,2OOO), 

/2(A,2) = 0(^,1100), 

/3(A,2) = 200000) 100100)}/^, 

/4(Ai,2) = {V'iA, 200000) + -0(^1, 100100) --0(^1, 110000)} /V3, 

MA,,2) = ^ 0«(F2,1OOO)V'^')(F2,1OOOOO)/V3, 
i/=i 

ME, 2) = 0^^)(E,11OO), 

/2(£;,2) = Ui'\E,200000)-4'\E,mi00)}/V2, 

fsiE, 2) = 200000) + 100100) - 2ipi^\E, 110000)} /V6, 

U{E,2) = 0(Ai, 1000)4') (S, 100000), 

ME, 2) = {0i')(F2, 1000)V'?)(F2, 100000) - 0^')(F2, mO)4'\F2, 100000)} /V2, 

/i(F2,2) = 0i')(F2,2OOO), 

/2(F2,2) = 0i')(F2,llOO), 

/3(F2,2) = {^^'^^2, 200000) -4^^^2,110000)}/^, 

/4(F2,2) = 4=^^(^2,110000), 

/5(F2,2) = 0«(F2, 1000)4') (E, 100000), 

/6(F2,2) = 0(Ai, 1000)4') (F2, 100000), 

MF2: 2) = {0? )(F2, 1000)4')(F2, 100000) + 0^')(F2, 1000)4') (F2, 100000)} /x/2, 

/i(Fi,2) = 0(')(Fi, 110000), 

/2(Fi,2) = 0i')(Fi, 1000)4') (E, 100000), 

MFi, 2) = {0i')(F2, 1000)4'^ (F2, 100000) - 0^')(F2, 1000)4^^(^2, 100000)} / A 

where the number 2 in the argument denotes v = 2. Those states belonging to the same 
irreducible representation will be mixed by the Hamiltonian. 

Directly calculating the Hamiltonian matrix H in those bases, we obtain a block matrix 
with 5x5 submatrix for 74i, 5x5 submatrix for £',7x7 submatrix for F2, and 3x3 
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submatrix for Fi: 



H{A,,2) = 



/ 2u:s - Qxs 





V 2V3X4 



2V3XI 

2ujs — 4xs + 4Ai 



-2A4 






Ci 
-^2/3x6 
4A4 



2^/3A4 \ 

-2A4 

/2j2>Xb 4A4 

C2 

C3 J 



H{E,2) 





/ 20;^- 


Axs - 2Ai 











4A4 \ 













Ci 




4A4 













-x^/a/S 


C4 

























C5 + 2A3 











l 


4A4 


4A4 








C3 ) 






( 2uJs — 6xs 


2Ai 








2A4 


2V2A4 






2Ai 


2us — 4X5 








2A4 


-2V2A4 

















2V2A4 




















C7 





4A4 












2V2A4 





C3 + 2A3 











2A4 


2A4 








C5 









2V2A4 


-2V2A4 





4A4 





C3 


} 



/ 2a;6-4x6-2A2 -2V2X4 \ 
H{Fi,2) = -2V2A4 C3 + 2A3 

V Cs J 

where Ci = 2ub — 5xb — 2A2 — 2A3, C2 = 2ujb — (14/3)xb — 2A2 + 2A3, C3 = Us — 2xs + uJb — 
2xfc-Ai-A2-A3, C4 = 2c^fe-(13/3)xfc-2A2 + 2A3, C5 = a;5-2a;5 + a;6-2a;6 + 3Ai-A2-A3, 
Ce = 2uJb - (16/3)a;b - 2A2, and C7 = 2ujb - 4:Xb - 2A2 - 2A3. 



IV. PURE VIBRATIONAL SPECTRA 



Now, we are going to fit tfie experimental data by our boson realization model. 
Methane is a typical molecule with symmetry. To our knowledge, there are 4 data 

for V = 1, 7 data for v = 2 and 8 data for v = 3. In our first scheme we fit those 19 data 
to determine the eight parameters (Table 1, the first scheme) with the root-mean-square 
energy deviation 11.61cm~^, where the standard deviation is calculated unweightedly: 

"1^ 19 

= E (^^(calc) - ^.(expt))' (4.1) 

Prom the results we come to three conclusions: 

i) Since Xb is relatively small, the bending oscillators are near harmonic ones. 

ii) The interbond interactions between bending vibrations are quite weak. 

iii) The interactions between stretching and bending vibrations are strong. 

From these conclusions, we propose our second scheme where the bending oscillators 
are harmonic {xb = 0) and there is no interaction between the bending vibrations (A2 = 
A3 = 0). The second scheme provides a five-parameter fit to the experimental data of 
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methane with the root-mean-square deviation 12.42cm ^ (see Table 1, the second scheme). 
Recall that Ref. [4] presented a seven-parameter fit with the root-mean-square deviation 

12.16cm-i 

For comparison, we list in Table 2 the 19 experimental data, the calculation results 
from the algebraic model [4], and our results in two schemes for the vibrational spectra 
{v < 3) of methane. 



Table I. Parameters in the Hamiltonian obtained by the least square fitting (cm ^) 



Scheme 


Stretching 


Bending 


Intcrac. 




ujs Xg Ai 


Uf, Xf, A2 A3 


A4 


rms 


First 


2986.74 77.96 34.55 


1508.37 -6.635 -5.96 -0.90 


-203.73 


11.61 


Second 


2986.24 76.55 33.60 


1525.85 


-201.65 


12.42 



Table II. Experimental data [8], algebraic calculation [4], and our calculation 
results for the vibration spectra {v < 3) of methane (cm~^). 



V = 1 


r 


F2 


E 


Ai 


F2 


Expt. [8] 


1310.0 


1533.0 


2916.5 


3019.4 


Calc. [4] 


1303.7 


1520.4 


2918.4 


3027.2 


Scheme 1 


1305.7 


1526.7 


2934.5 


3019.0 


Sclieiiie 2 


1307.8 


1525.9 


2933.9 


3017.6 



V = 2 and F2 


Expt. [8] 


2614.0 


2830.4 


4223.0 


4319.0 


4549.0 


5861.0 


6004.7 


Calc. [4] 


2610.5 


2841.5 


4222.0 


4330.9 


4547.7 


5856.7 


6014.5 


Scheme 1 


2610.1 


2840.1 


4226.7 


4308.4 


4546.9 


5855.8 


6011.2 


Scheme 2 


2614.3 


2833.6 


4228.8 


4309.5 


4543.5 


5855.3 


6008.6 



V — 3 and F2 


Expt. [8] 


4123.0 


5775.0 


5861.0 


7514.0 


8604.0 


8807.0 


8900.0 


9045.0 


Calc. [4] 


4123.9 


5759.9 


5868.7 


7534.9 


8603.0 


8794.1 


8910.0 


9034.5 


Scheme 1 


4136.9 


5759.3 


5858.4 


7513.8 


8601.9 


8805.3 


8915.9 


9035.7 


Scheme 2 


4140.1 


5754.6 


5851.2 


7534.4 


8603.3 


8804.1 


8913.0 


9031.4 



In terms of the eight parameters in the first scheme or the five parameters in the 
second scheme, it is straightforward to calculate the rest of the vibrational spectra for 
methane. After removing the spurious states, for v = 2, there are 5 states with Ai, 5 sets 
of states with E, 7 sets of states with F2, and 3 sets of states with Fi. For v = 3, there 
are 13 states with Ai, 4 states with A2, 14 sets of states with E, 25 sets of states with F2, 
and 15 sets of states with Fi. Except for the states with F2 and v — 2 that were listed in 
Table II, the rest of calculation results in the first scheme are listed as follows. 

v = 2, Ai 2614.3 3057.9 4297.7 5807.6 5974.4 
v^2, E 2616.5 3055.6 4326.3 4461.2 6038.2 
v^2, Fi 2832.4 4324.8 4545.7 
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V. CONCLUSIONS 

In this paper we describe ten coupled one- dimensional anharmonic oscillators of a 
tetrahedral molecule by ten sets of bosonic creation and annihilation operators. The ten 
oscillators are divided into two classes: stretching and bending oscillators. The energy 
levels of those oscillators are described by four parameters under the assumption of the 
Morse potential for stretching vibration and the Poschl- Teller potential for the bending 
vibrations: uUs, Xs, lvi,, and Xi,. The intcrbond interactions and the interactions between 
stretching and bending vibrations are supposed to be invariant and to preserve the 
total number v of vibrational quanta so that 4 parameters Aj, 1 < i < 4, have to be 
introduced. 

In the first scheme of the boson realization model with eight parameters we fit the 
19 experimental vibrational data for methane, and obtain the root-mean-square energy 
deviation to be 11.61 cm~^. Prom the obtained parameters, we see that the interactions 
between the bending vibrations are weak, the interactions between the stretching and 
bending vibrations are strong, and the bending oscillators are quite near harmonic ones. 
These conclusions are different from the previous model [4]. From these conclusions 
we proposed another five-parameter fit in the second scheme with the root-mean-square 
energy deviation 12.42cm~^. To our knowledge, it may be the model with the least 
parameters that well fits the experimental vibration spectra {v < 3) of methane. 

The interaction between vibrational and rotational motions plays an important role in 
describing the abundant experimental data of vibrorotational energy spectra of a tetra- 
hedral molecule. We will study it by the boson realization model elsewhere. 
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